26.2 SCOPE vs STOP

Fit LME models for Prolactin, HPL and GH2.
Fit GEE models for Insulin and HOMA-IR.

  • Outcomes: Prolactin, hPL, GH2, Insulin, HOMA-IR, Glucose (all outcomes log-transformed)
  • Independent variables: Study, GA at sampling
  • Model fitted: LME or GEE model
          log(hormones) ~ Study + GA at sampling + Study:GA at sampling
          LME only: Random intercepts for Study ID and Plate
          GEE only: Independence correlation structure

26.2.1 Global tests

Type III ANOVA on main and interaction effects, and LM model summaries.

PRL (ng/ml)
Chisq Df Pr(>Chisq)
(Intercept) 8.706512 1 0.003 **
Study 2.244395 1 0.1 ns
GA_samp 30.790407 1 <0.0001 ***
Study:GA_samp 2.221336 1 0.1 ns
hPL (μg/ml)
Chisq Df Pr(>Chisq)
(Intercept) 46.94900 1 <0.0001 ***
Study 21.66283 1 <0.0001 ***
GA_samp 94.21469 1 <0.0001 ***
Study:GA_samp 25.95307 1 <0.0001 ***
GH2 (ng/ml)
Chisq Df Pr(>Chisq)
(Intercept) 12.13572 1 0.0005 ***
Study 16.86104 1 <0.0001 ***
GA_samp 53.71079 1 <0.0001 ***
Study:GA_samp 17.66908 1 <0.0001 ***
Insulin (uU/ml)
model term df1 df2 F.ratio Chisq p.value
1 Study 1 Inf 6.440 6.440 0.01 *
3 GA_samp 1 Inf 3.036 3.036 0.08 ns
2 Study:GA_samp 1 Inf 4.677 4.677 0.03 *
HOMA-IR
model term df1 df2 F.ratio Chisq p.value
1 Study 1 Inf 7.472 7.472 0.006 **
3 GA_samp 1 Inf 3.420 3.420 0.06 ns
2 Study:GA_samp 1 Inf 4.874 4.874 0.03 *
Glucose (mmol/L)
model term df1 df2 F.ratio Chisq p.value
1 Study 1 Inf 0.232 0.232 0.6 ns
3 GA_samp 1 Inf 0.008 0.008 0.9 ns
2 Study:GA_samp 1 Inf 0.111 0.111 0.7 ns

26.2.2 Model summary

Ratio of Geometric means (95% CI)

PRL (ng/ml)
Var Effect P
(Intercept) 1.72 (0.58, 2.87) 0.003
StudySTOP 0.91 (-0.28, 2.10) 0.13
GA_samp 0.21 (0.13, 0.28) <0.0001
StudySTOP:GA_samp -0.06 (-0.14, 0.02) 0.14
hPL (μg/ml)
Var Effect P
(Intercept) 2.14 (1.52, 2.76) <0.0001
StudySTOP -1.51 (-2.15, -0.87) <0.0001
GA_samp 0.20 (0.16, 0.24) <0.0001
StudySTOP:GA_samp 0.11 (0.07, 0.15) <0.0001
GH2 (ng/ml)
Var Effect P
(Intercept) -1.33 (-2.10, -0.56) 0.0005
StudySTOP -1.64 (-2.45, -0.83) <0.0001
GA_samp 0.18 (0.13, 0.23) <0.0001
StudySTOP:GA_samp 0.11 (0.06, 0.17) <0.0001
Insulin (uU/ml)
Var Effect P
(Intercept) 4.31 (2.75, 5.87) <0.0001
StudySTOP -2.10 (-3.73, -0.48) 0.01
GA_samp -0.11 (-0.21, -0.01) 0.04
StudySTOP:GA_samp 0.12 (0.01, 0.23) 0.03
HOMA-IR
Var Effect P
(Intercept) 3.01 (1.39, 4.64) 0.0003
StudySTOP -2.36 (-4.05, -0.67) 0.006
GA_samp -0.12 (-0.22, -0.01) 0.03
StudySTOP:GA_samp 0.13 (0.01, 0.24) 0.03
Glucose (mmol/L)
Var Effect P
(Intercept) 1.67 (1.34, 1.99) <0.0001
StudySTOP -0.08 (-0.42, 0.26) 0.63
GA_samp 0.00 (-0.02, 0.02) 0.90
StudySTOP:GA_samp -0.00 (-0.03, 0.02) 0.74

26.2.3 Diagnostics

Model diagnositcs.

Diagnositcs for PRL (ng/ml)

Diagnositcs for hPL (μg/ml)

Diagnositcs for GH2 (ng/ml)

Diagnositcs for Insulin (uU/ml)

Diagnositcs for HOMA-IR

Diagnositcs for Glucose (mmol/L)

26.2.4 Estimated means

Estimated marginal means (or adjusted means) adjusted for GA at sampling.

PRL (ng/ml)
Study GA_samp response SE df lower.CL upper.CL
SCOPE 10 44.67877 9.765118 775.60673 29.09173 68.61720
STOP 10 60.90187 4.717695 79.27188 52.20009 71.05425
SCOPE 16 155.26985 12.165875 82.42309 132.86150 181.45757
STOP 16 147.84441 14.552006 197.53829 121.76053 179.51605
hPL (μg/ml)
Study GA_samp response SE df lower.CL upper.CL
SCOPE 10 60.06118 6.777997 1013.3294 48.13041 74.94940
STOP 10 39.61791 1.055110 108.7181 37.58090 41.76535
SCOPE 16 194.04839 5.259046 113.7242 183.90473 204.75155
STOP 16 247.31102 10.211412 486.6536 228.03943 268.21126
GH2 (ng/ml)
Study GA_samp response SE df lower.CL upper.CL
SCOPE 10 1.600762 0.2192977 433.4993 1.222895 2.095387
STOP 10 0.948957 0.0264086 216.5576 0.898308 1.002462
SCOPE 16 4.703807 0.1314798 197.9604 4.451542 4.970366
STOP 16 5.462550 0.2814179 922.9232 4.937259 6.043728
Insulin (uU/ml)
Study GA_samp response SE df asymp.LCL asymp.UCL
SCOPE 10 25.31505 7.1937998 Inf 14.504141 44.18404
STOP 10 10.17569 0.4433112 Inf 9.342878 11.08274
SCOPE 16 13.27781 0.7140309 Inf 11.949565 14.75370
STOP 16 10.90728 1.0191172 Inf 9.082069 13.09930
HOMA-IR
Study GA_samp response SE df asymp.LCL asymp.UCL
SCOPE 10 6.371755 1.8874508 Inf 3.565455 11.386841
STOP 10 2.137866 0.0979181 Inf 1.954312 2.338660
SCOPE 16 3.173325 0.1800916 Inf 2.839275 3.546678
STOP 16 2.273611 0.2215450 Inf 1.878337 2.752067
Glucose (mmol/L)
Study GA_samp response SE df asymp.LCL asymp.UCL
SCOPE 10 5.361853 0.3217762 Inf 4.766861 6.031111
STOP 10 4.750210 0.0569715 Inf 4.639850 4.863194
SCOPE 16 5.406086 0.0560985 Inf 5.297246 5.517163
STOP 16 4.682279 0.0697641 Inf 4.547522 4.821030

26.2.5 Comparisons

Hormone Contrast GA_samp Ratio of Geometric means (95% CI) P-value
PRL (ng/ml) STOP / SCOPE 10 1.36 (0.90, 2.07) 0.1
PRL (ng/ml) STOP / SCOPE 16 0.95 (0.81, 1.11) 0.5
hPL (μg/ml) STOP / SCOPE 10 0.66 (0.53, 0.82) 0.0003
hPL (μg/ml) STOP / SCOPE 16 1.27 (1.18, 1.38) <0.0001
GH2 (ng/ml) STOP / SCOPE 10 0.59 (0.45, 0.78) 0.0002
GH2 (ng/ml) STOP / SCOPE 16 1.16 (1.04, 1.29) 0.007
Insulin (uU/ml) STOP / SCOPE 10 0.40 (0.23, 0.71) 0.002
Insulin (uU/ml) STOP / SCOPE 16 0.82 (0.67, 1.01) 0.07
HOMA-IR STOP / SCOPE 10 0.34 (0.19, 0.60) 0.0003
HOMA-IR STOP / SCOPE 16 0.72 (0.57, 0.89) 0.003
Glucose (mmol/L) STOP / SCOPE 10 0.89 (0.79, 1.00) 0.05
Glucose (mmol/L) STOP / SCOPE 16 0.87 (0.84, 0.90) <0.0001

26.2.6 Barplots by Study

Barplots in this section may be misleading due to involvement of Study by GA interaction.

26.2.7 Plots for GA